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Marginally unstable Holmboe modes for smooth density and velocity profiles are studied. For a 
large family of flows and stratification that exhibit Holmboe instability, we show that the modes 
with phase velocity equal to the maximum or the minimum velocity of the shear are marginally 
unstable. This allows us to determine the critical value of the control parameter R (expressing 
the ratio of the velocity variation length scale to the density variation length scale) that Holmboe 
instability appears R cr u = 2. We then examine systems for which the parameter R is very close to 
this critical value R cr u- For this case we derive an analytical expression for the dispersion relation 
of the complex phase speed c(fc) in the unstable region. The growth rate and the width of the 
region of unstable wave numbers has a very strong (exponential) dependence on the deviation of R 
from the critical value. Two specific examples are examined and the implications of the results are 
discussed. 

PACS numbers: 



I. INTRODUCTION 

Holmboe instability in stratified shear flows appears 
in a variety of physical contexts such as in astrophysics, 
the Earth's atmosphere and oceanography 0, 0,11, M 
EI IE El Although the typical growth rate is smaller 
than the one of Kelvin-Hclmholtz instability it is present 
for arbitrarily large values of the global Richardson num- 
ber making Holmboe instability a good candidate for the 
generation of turbulence and mixing in many physical 
scenarios. 

What distinguishes Holmboe from the Kelvin- 
Helmholtz instability is that unlike the later instability 
the Holmboe unstable modes have non-zero phase veloc- 
ity that depends on the wavenumber (i.e. traveling dis- 
persive modes). It was first identified by Holmboe 
in a simplified model of a continuous piece-wise linear 
velocity profile and a step-function density profile. Sev- 
eral authors have expanded Holmboe's theoretical work 
[Til IrH Il3l Il4j by considering different stratification and 
velocity profiles that do not include the simplifying sym- 
metries Holmboe used in his model. Hazel and more 
recently Smyth and Peltier 0] and Alexakis El have 
shown that Holmboe's results hold for smooth density 
and velocity profiles as long as the length scale of the 
density variation is sufficiently smaller than the length 
scale of the velocity variation. Furthermore, effects of 
viscosity and diffusivity El 0, non-linear evolution 
[U HI Hi E3 and mixing properties [H of the Holmboe 
instability have also been investigated. The predictions of 
Holmboe have also been tested experimentally. Browand 
and Winant [2{| first performed shear flow experiments 
in a stratified environment under conditions for which 
Holmboe's instabilities are present. Their investigation 
has been extended further by more recent experiments 

El El El El IH m IH E3 . 

Although the understanding of Holmboe instability has 
progressed a lot since the time of Holmboe, there are still 



basic theoretical questions that still remain unanswered, 
even in the linear theory. Most of the work for the linear 
stage of the instability has been based on the Taylor- 
Goldstein equation (see H^|), which describes linear nor- 
mal modes of a parallel shear flow in a stratified, inviscid, 
non-diffusive, Boussinesq fluid: 



dy 2 



U" 



J(y) 



U-c (U-c) 2 



= 0, 



(1) 



where 4>(y) is the complex amplitude of the stream func- 
tion for a normal mode with real wavenumber k. c is 
the complex phase velocity. Im{c} > implies instabil- 
ity with growth rate given by C — klm{c}. U(y) is the 
unperturbed velocity in the x direction. J(y) — —gp'/p 
is the squared Brunt- Vaisala frequency where p is the 
unperturbed density stratification and g is the accel- 
eration of gravity. Prime on the unperturbed quanti- 
ties indicates differentiation with respect to y. Equa- 
tion Q together with the boundary conditions <fi — > 
for y — ■* ±oo, forms an eigenvalue problem for the com- 
plex eigenvalue c. Here we just note a few known re- 
sults for the Taylor-Goldstein equation. If c is real 
and in the range of U there is a height y c , at which 
U{y c ) = c. At this height y Cl called the critical height, 
equation Q has a regular singular point. For some con- 
ditions unstable modes exist with the real part of the 
phase velocity within the range of U . The phase speed 
of these modes satisfies Howard's semi-circle theorem 
|c - l/2(sup{?7} + inf{J7})| < 1/2 j sup{[/} - inf{!7}|. 
If these unstable modes exist the Miles-Howard theo- 
rem [34j guarantees that somewhere in the flow the local 
Richardson number defined by: 



My) 



J(y) 
[U'(y)Y 



(2) 



must be smaller than 1/4. 

A typical example used in many studies assumes 
a velocity profile given by U(y) = tanh(y) and the 
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FIG. 1: Stability diagram for the Hazel model for R — 3. 
Panel (a) shows the Kelvin-Helmholtz unstable regions with 
dark gray, and Holmboe unstable regions with light gray. 
Panel (b) shows the same instability diagram for larger values 
of Jo. 



squared Brunt- Vaisala frequency being given by J = 
Jo cosh(Ry)~ 2 . Where Jo is the global Richardson num- 
ber usually denned as Jo = Ri(0). The case of R = 1 was 
examined by Miles |3fij analytically and numerically by 
Hazel ^f|. It exhibits only Kelvin-Helmholtz instability 
for the wavenumbers that satisfy k(l — k) > Jq. Hazel 
|15| also examined numerically the case R = 5 where it 
was shown that along with the Kelvin-Helmholtz unsta- 
ble region there is also a stripe (in a Jo — k diagram) of 
unstable Holmboe modes. As an example the instability 
region for R = 3 is shown in figure ^ Hazel observed 
that if R > 2 there is always a height at which the local 
Richardson number is smaller than 1/4. Based on this 
observation Hazel conjectured that R = 2 is the critical 
value of R above which the Holmboe instability appears. 
Later careful numerical examination by Smyth jlq j found 
unstable Holmboe modes to appears only for values of R 
larger than R > 2.4. More recently Alexakis ^3 showed 
that the instability can be found for smaller values of R 
up to R — 2.2 making the conjecture by Hazel still plausi- 
ble. More specifically Alexakis 01 showed (numerically) 
for the examined shear and density profiles that the left 
instability boundary of the Holmboe instability region 
(see figure^,) is composed of marginally unstable modes 
with phase velocity equal to the maximum or the min- 
imum of the shear velocity. Such a condition has been 



known to hold for smooth velocity and discontinuous den- 
sity profiles [36l l37l l38| . Finding these marginally unsta- 
ble modes corresponds to solving a Schrondiger problem 
for a particle in a potential well: 



where 



dy 2 



v c {y) 



[k 2 + V c (y)] = 0, 



u" J(y) 

U~c {U-cf 



(3) 



(4) 



and c is taken to be c = U, 



max / min 



the maximum or min- 



imum velocity of the shear layer. The right boundary of 
the Holmboe unstable region on the other hand (see fig- 
ure^) is composed of singular modes with phase velocity 
within the range of the shear velocity. These modes can 
be determined by imposing the condition that the solu- 
tion close to the critical height can be expanded in terms 
of only one of the two corresponding Frobenius solutions. 
Furthermore it was shown that for sufficiently large Jo 
more than one instability stripe exists. These new insta- 
bility stripes are related with the higher internal gravity 
modes of the unforced system. Figure \T]p shows the in- 
stability region for R = 3 and Jo up to 80. We note that 
different unstable Holmboe modes have been found ex- 
perimentally in |32| that were then interpreted in terms 
of the multi-layer model of [T^ . 

The understanding however of the linear part of Holm- 
boe instability for smooth shear and density profiles still 
remains conjectural and most of the results are based on 
numerical calculations and do not therefore constitute 
proofs. We try to address some of these issues in the 
present work. In the next section we prove for a gen- 
eral class of velocity profiles that the modes that have 
phase velocity equal to the maximum/minimum velocity 
of the shear are marginally unstable. In section II I II we 
examine the case for which the parameter R is slightly 
larger than it's critical value, and the dispersion rela- 
tion inside the instability region is derived based on an 
asymptotic expansion. In section lTVI we test these results 
for specific shear and density profiles. A summary of the 
results and final conclusions are in the last section. To 
help the reader, a table [I] with the definitions of all the 
basic symbols used in this paper has been added. 



II. MARGINAL WAVENUMBER 

In this section we examine modes with phase veloc- 
ity equal to the maximum/minimum velocity of the flow, 
and show under what conditions these modes constitute 
a stability boundary. We start by considering an infi- 
nite shear layer specified by the monotonic velocity pro- 
file U(y) that has the asymptotic values U(±oo) = L/±oo- 
Since the system is Galilean invariant with no loss of gen- 
erality we can set U +OQ = — U-oo = U^. More precisely 
we will assume that the asymptotic behavior of U(y) for 
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Symbol 


Definition 


c 


Complex phase velocity 


c 


= fclm{c} the growth rate 


k 


Wavenumber 


<Kv) 


Stream function 


k 


Wavenumber for which c(fco) = Uoo 


4>o 


Stream function for the c(fco) = Uoo mode 


Vc 


Location of the critical layer U(y c ) = c 


Jo 


Global Richardson Number 


Uoo, U*,a 


Coefficients that appear in the large y 




behavior of U(y) ~ U x - U*e~ ay 


P,J* 


Coefficients that appear in the large y 




behavior of J(y) — J*e~ l3y 


<fioc, A 


Coefficients that appear in the large y 




behavior of 4>o{y) — ifiooe~ Xv 


Joe, cr, R, q 


J*/{U*af, U*/Uoe, a/0, k/a 


5 


= R~2 


J 


=Joo(e Cl /a) s 




local Richardson number at y = y c 




= 1/2-^1/4-7 


A 


= \Jq 2 + 1 — J 


F(a, b, d, s) 


Hypergeometric Function 



TABLE I: Table of used symbols 

y — > +00 is going to be given by 

U[y) ~ C/oo - U*e- ay . 

We further assume that the layer is stably stratified with 
J(y) > having asymptotic behavior J(y) ~ J*e~P y for 
y — ► +00. In what follows we are going to concentrate 
only on the modes with phase velocity close to c = Uoo', 
the results can easily be reproduced for the c = U-oc 
modes by following the same arguments. To simplify the 
problem we will non-dimentionalize the equations using 
the maximum velocity Uoo and the length-scale a -1 (es- 
sentially setting Uoo = 1 and a = 1). The resulting non- 
dimensional control parameters for our system are the 
asymptotic Richardson number: Joo = J*/(i7*a) 2 , the ra- 
tio of the two velocities:*! = U*/Uoo, and the ratio of the 
two length scales: R = (3 /a. The wavenumber becomes 
q = k/a and c is measured in units of Uoo- 

Finally we assume that a solution 4>o(y) of the 
Schrondiger problem described in equation [3] for the 
wavenumber qo — ko/a exists. Clearly if R > 2 and c = 1 
the asymptotic behavior of V c for large y is V c (y) ~ 1 and 
(j>o(y) behaves as <j> ~ (^ 00 e _Ay with A = y/q 2 + 1. If how- 
ever we have R = 2 then V c (y) ~ 1 — Joo for y — * 00 and 



a = vVTT 

both cases 



Joo- For abbreviation we will denote for 



A= V q 2 + 1 - J 



where J (that will be defined precisely later on) takes 
the values J ~ when R > 2 and J = Joo when 
R = 2. As discussed in 01 j no solution exists that satis- 
fies the boundary conditions for the Schrondiger problem 




FIG. 2: The potential V c for the Hazel model with R = 3 
and Jo = 1 for three values of c, c = 1 (solid line), c = 1 + e 
(dashed line), c = 1 — e (dotted line), where < £ < 1. For 



c = 1 V c asymptotes to the value V c 
drastically changes when c / 1. 



a . This behavior 



described in equation [3] if R < 2 since it corresponds in 
finding bounded eigenstates in an unbounded potential 
well. 

Our aim in this section is to find how c changes from 
the value 1 as we increase q from the value qo — fco /ex. We 
proceed by carrying out a regular asymptotic expansion 
by letting q = go + eqi and c = 1 — eci + . . . with e«l 
and c\ in general complex. However as we deviate from 
the c = 1 solution the behavior of the potential V c (y) 
drastically changes (0(1) change) in the large y region 
and only slightly (linearly with respect to the change in 
c) for y ~ 0(1). Figure [3 illustrates this change. This 
implies that two different expansions are needed, one for 
y being of 0(1) and one for large y. 



A. Local solution: y — 0(1) 

We begin with the local solution and expand <f> as cj) — 
b + £0i + . . . . For y = 0(1) at first order we obtain: 

[ql + V (y)] 0i = [2gig - ctVtiy)] O (5) 



dy 2 



where Vb = V c is given in equation 0] for c = 1 and 



Vi(y) = 



d_ 
~dc 



U" 



u" J(y) 

U-c (U-c) 2 
2J(y) 



c=l 



(Z7- 1) 2 (C/-1) 3 ' 

The solution of this inhomogeneous equation can be 
found using the Wronskian to obtain: 



4>i = My) 



dy' (6) 



4 



Where we have used as normalization condition 0(0) = 
0o (0). Clearly this solution satisfies the boundary con- 
dition for y — > -co. For y — > +co, by performing the 
integrations we obtain: 



2q 1 q I 1 - ah 



-t Xy + 0(e- Xy ) 



(7) 



where I x = /_ ™ 0gdy > and J 2 = /_ ~ V^dy. So the 
large y behavior of based on the local solution is given 
by: 



2q 1 q I 1 - al 



-Xy , ^yiyu^i 



(8) 



B. Far away solution: y = C(ln[l/e]) 

To capture the large y behavior we need to make the 
change of variables y — y — y c where y c = — hi(eci/er) is 
the location of the singularity determined by U(y c ) = c. 
(Note that for c\ complex y does not coincide with the 
real y axis.) The Taylor-Goldstein equation^then reads: 



d 2 ^ 
dy 2 



e» - 1 



(ev - l) 2 



= (9) 



where 6 = R — 2 and only the leading terms have been 
kept. Introducing the variable s — e~ y we obtain: 



the solution of which is the Hypergeometric function 
h(s) = F(a, 6, d; s) with: 



a= (fJL + q) + Jq 2 + 1-J, 



(m + q) - yq 



i- J, 



and 



(2-7+1) 



Note that q + /i — a = —A and q + /i — b = +A. Some 
basic properties of the Hypergeometric function are given 
in appendix A, here we give just some of the resulting 
asymptotic behavior of <f>: 



lim ( 

lim < 

s— >+oo 



o9 — „-qy 



T(d)T(b-a) 
T(b)T(d-a) [ ' 
T(d)T(a - b) 
T(a)T(d-b) 



s«(-s) 



fj,—b 



(13) 



(14) 



Returning to the y variable and up to a normalization 
factor A we have for y <C y c that 



.4 



-a, , ( )2X T(a)T(d-b)T(-2X) Xv 
+ { l) T(b)T(d-a)T(2\) 



(15) 



ds 2 



ds 



1 



Js 2 g* 
(1-s) 2 



= 0, (10) 



where J = ■ (ecx/cr) gives the Richardson number at 
the critical height. Note that if R = 2 (i.e. S = 0) then 
J = Joo = 0(1). If however R > 2 then J < 1 and the 
term in the brackets proportional to J is small and can 
be neglected everywhere except close to the singularity 
s — 1 . To deal with this small singular term we can write 
for s close to 1: s 5 ~ 1 — 5(1 — s) + . . . and keep the leading 
term. That way the principal term inside the brackets is 
always kept for all values of s and our solution will be 
correct to first order for all values of R > 2 by solving: 



C. Matching 

Matching the exponentially decreasing terms of the lo- 
cal and the far-away solution we obtain: A = tpooi an d 
from the exponentially increasing terms we have: 

2<Mo/i -c x h . ,2X T(a)T(d - b)T(-2X) 



•■fo 



T(b)T(d-a)T(2\) 



(16) 

We can solve the equation above iteratively by letting 
eci = ec[ + e 2X c' 2 + . . . . To first order we obtain: 



2qiq h 



(17) 



ds 2 



d 
ds 



Js 2 



1-s (1-s) 



= 0. (11) 



To deal with the singularities at s = and s = 1 we 
make the substitution = s q (l — s) fl h(s) with /.i = 1/2 — 

1/4 — J. We then obtain: 



S (l 



S) d^ h 



[(2q+l) - (2 f i + 2q 



[l-/i-( ? + l)]ft = 0, 



(12) 



which gives the first correction to the phase speed and 
determines if the real part of phase speed is increasing or 
decreasing with the wavenumber. If for example I2 > 
(which will be the case in the examples that follow) then 
c is decreasing with q and the correction c[ is positive for 
positive qi and negative for negative q\. The opposite 
holds if I2 < 0. From now on we will assume that I2 > 
which is the physically relevant case (Re{c(fc)} being a 
decreasing function of k) if however there is a velocity 
and density profile such that I2 < the same results will 
hold but for the opposite direction in q (i.e. wavenumbers 
smaller than qo will be unstable and wavenumbers larger 
than qo will be stable). The correction however is real, 



■5 



and contains no information about the growth rate. At 
the next order we have 



T(a)T(d-b)n-2X) 
Voo( Cl ) r(6)r(d _ a)r(2A) ■ (18) 



This correction is much smaller but contains the first or- 
der correction of the imaginary part of c. The dispersion 
relation of c for q close to go can then be written in terms 
of g as: 



c= 1- 



2q a L 



ion, \ 



h 



-2koh 
h 



(q - go) 



2\ 



T{a)T{d-b)T(~2X) 
T(b)T(d-a)T(2X) 



(19) 



Special care is needed to interpret the term (— c[) for 
c[ is given by equation El When qi < 0, c[ is nega- 
tive and the term (— c' 1 ) 2A is real, this corresponds to the 
case that c becomes larger than the shear velocity and no 
critical layer is formed. The Howard semi-circle theorem 
then guaranties stability. This proves that wavenumbers 
slightly smaller than q are stable. When q\ > 0, c[ 
is positive and (— c[) 2X becomes a complex number that 
can take different values depending on whether the mi- 
nus sign is interpreted as e l7r or e~ l7V . The choice de- 
pends on the location of the singularity on the complex 
plane when we integrate the Taylor Goldstein equation 
^ If Im{ci} > then (— ci) 2A should be interpreted 
as | ci | 2 A e~ l2A7r because the integration is going over the 
singularity. If Im{ci} < then (— ci) 2A should be inter- 
preted as |ci| 2A e +l2A7r because the integration is going 
under the singularity. Here we arrive at an important 
point in our derivation: the sign of the imaginary part of 
c based on equation 1191 depends on the original assump- 
tion about the sign of Im{c} when we integrate across 
the singularity. Thus, in order for the matching to be 
successful we need to verify that the original assumption 
about the sign of Im{ci} is consistent with the final re- 
sult. If we assume that Im{ci} > then fromEl we have 
that 



< Im{ci} = sin(2A7r)|ci 



2A ^ 2 0O r( Q )r(d-b)r(-2A) 

1 I 2 T(b)T(d- a)T{2X) 



where Im{(c' 1 ) 2A } is written as — sin(2A7r)|ci | 2A as pre- 
viously discussed. The matching is successful only if the 
sign of the right hand side (r.h.s.) of eauationl2()lis posi- 
tive as originally assumed and only then is the dispersion 
relation E| valid. (We arrive at the same condition if we 
initially assume that Im{c} < 0). It is shown in appendix 
B that for R > 2 (i.e. J ~ 0) the r.h.s. of|2U|is always 
positive and the matching is successful. For the special 
case however that R = 2 (i.e. J = J^) the matching is 
not always successful because the product T(b)T(d — a) 
that appears in equation 1201 can change sign depending 
on the value of J. In particular it is shown in the ap- 
pendix that if J > 2q/(2q + l) 2 the r.h.s. of equation 



1201 is negative and thus we end up with a contradiction. 
Therefore in the R = 2 case we have shown instability 
only if 



J<2g/(2g+l) 2 . 



(21) 



The unsuccessful matching when the condition is not 
satisfied suggests that the modes with q > qo are part of 
the continuous spectrum of the Taylor Goldstein equa- 
tion. This kind of modes have a discontinuity of the first 
derivative of <f> at the critical lay er and have been studied 
before in the literature [3^, |4£| . We need to emphasize 
here that the lack of instability at this order does not 
imply stability. Non-zero growth rate of smaller order 
can still exist and therefore the above result should be 
interpreted only as a sufficient condition for instability. 

To summarize we have shown that if R > 2 the modes 
with phase velocity equal to the maximum phase velocity 
of the shear (when they exist, and for density and velocity 
profiles that satisfy the conditions stated at the beginning 
of this section) are marginally unstable: wavenumbers 
with q < qo are stable and wavenumbers with q > qo are 
unstable. If R = 2 these modes are marginally unstable 
only if the condition is further satisfied and stable (to 
the examined order) otherwise. 



III. MARGINAL R 

In the last sectioned we showed marginal instability 
when the wavenumber q is varied from the critical value 
go- However the wavenumber is not a control parameter 
in a system. We would like therefore to examine a system 
for which one of the control parameters (Jo or R) is close 
to the critical value for which the instability begins. Since 
Holmboe instability is present for arbitrary large values 
of Jo the only other control parameter left is R. We con- 
sider therefore a case for which R = 2 + 5 with < 6 <C 1 
and the c = 1 solution (</>o,go) is known with go such 
that Joo > 2g /(2g + l) 2 so that the R = 2 case gives 
no instability at the examined order. We make a small 
variation in g = go + egi and c = 1 — ec\ with the exact 
relation between 8 and e still undetermined. At this stage 
we assume that e is sufficiently smaller than 8 so that the 
procedure in the previous section is still valid and then 
gradually increase its value until the approximations in 
the previous section start to fail. As we increase the value 
of e the most sensitive term (in e) that will be affected 
first, is the term proportional to J = J^ • (ec\/cr) S in 
the equation II II for which e is raised to the smallest ap- 
pearing power. Note that if e <C exp[— 1/8} then Jcl 
and the results of the previous section are still valid. If 
however e ~ 0(exp[— 1/8]) then J ~ 0(1). Following the 
same steps as in the previous section we end up in the 
dispersion relation given by equation 1191 but like in the 
R = 2 case J cannot be treated as a small parameter. 

The difference from the 8 = 0(1) case will therefore 
appear when we try to determine the sign of the r.h.s. 
of equation [201 To have successful matching we need to 
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satisfy the condition [2] Since J is finite the condition 
J < 2q/(2q + l) 2 that also appears in the R = 2 case 
could be violated. To capture the whole unstable region 
we define e such that Joo^ 5 = 2q , o/(l + 2g ) 2 or 



2q /J a 



L(l + 2g ) 2 



1/8 



< 1. 



(22) 



Note that the term inside the brackets is always smaller 
than one. For such a choice the condition [21 for instabil- 
ity reads 



J = J> 



/eci 
V a 



2<?o 



(1 + 2q ) 2 



[f + <51n( Cl /a) +0(6 2 )} 
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k 
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FIG. 3: The c = 1 solutions for the Hazel model Jo(fc) for R = 
2 (solid line), the condition from equation l21I Jn < 4fc/(fc+l) 2 
(dashed line). 



< 



2go 

(i + 2 qo y 



O(e) 



(23) 



IV. EXAMPLES 



or ci/er < 1. Already at this stage it can be seen that we 
have instability only if c\ = 2q\qoI\j I2 < o and therefore 
the instability is confined in the region of wave numbers 



qo < q < q + Aq 



(24) 



where Aq = eal2/2q Ii. Therefore, the second insta- 
bility boundary for the Holmboc instability is given by 
q + Aq. To get the full dispersion relation in this asymp- 
totic limit we need to expand in terms of S the product 
T(d — a)T(b) that appears in equation EH since this is 
the term that can change sign depending on the value of 
J. This is done in Appendix B and the resulting growth 
rate inside the instability region to the first non-zero or- 
der becomes: 



/- t r l xr< \ i2A 1 f 2(<7 - <7o)<?o^i 

C = <7oIm{c} = -ddqolqo - q\ In 1 



(25) 

where C\ > is an 0( I) quantity and is given in equation 
IBII The maximum of the growth rate is obtained for 
q — qo = €e~ x t 2X I20 / {2qall) with the growth rate being 
given by: 



max [£] = Se 



2Xe 



he 



2X 



(26) 



Therefore the growth rate scales like e 2A and the width 
of the instability region scales like Aq ~ e. In terms of 5 
these relations are given by £ ~ Se~ 2Xl ^ s and Aq ~ e -7 / 5 
where 7 is a positive constant. This very strong depen- 
dence with 5 suggests that both £ and Aq decrease very 
rapidly as S becomes smaller. This can explain the diffi- 
culty numerical codes have, when attempting to calculate 
growth rate for values of R very close to R = 2. 



In the previous sections we showed some general re- 
sults for the Holmboe unstable modes. In this section we 
examine some specific examples often used in the litera- 
ture to model Holmboe's instability. We begin with the 
model introduced by Hazel that we briefly mentioned in 
the introduction. The model assumes a velocity profile 
given by U (y) = tanh(y) and a density stratification de- 
termined by —gp'/p = Jo cosh -2 (Ry). Based on the def- 
initions given in section [H] we have that a = 2, j3 = 2R, 
U* = 2 and J* = 4 Jo. The resulting non-dimensional 
quantities are = Jo/4, a — 2, q — k/2 and R has 
the same meaning. This model satisfies all the condi- 
tions that are stated in section [H] therefore for R > 2 the 
modes with c(fc) = ±1 are marginally unstable. Further- 
more for the case R = 2 we will have instability only if 
the conditionl21lis satisfied, or in the units of this exam- 
ple if Jo/4 < k/(k + l) 2 . A simple numerical integration 
shows that this is not the case for this profile (see figure 
0. Therefore, the R — 2 is stable (to the examined or- 
der) and is the critical value beyond which the Holmboe 
instability begins. The imaginary part of c(k) for this 
profile for the case that R — 2.1 and Jo = 1.2 is shown if 
figure 01 where the numerical result is compared with the 
asymptotic expansion of equation 1251 Although S = 0.1 
is not very small there is satisfactory agreement (a 20% 
difference) between the asymptotic and the numerical re- 
sult. It is worth mentioning that it is very hard to find 
a range of values of S that both the asymptotic result 
is valid and Im{c} is large enough to be captured by a 
numerical code. Note that by decreasing the value of S 
from 0.1 to 0.05 has resulted in a drop of Im{c} by three 
orders of magnitude. 

A second family of flows we will examine assumes a 
velocity profile given by U(y) — tanh(y) as in the Hazel 
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FIG. 4: The Im{c} for the Hazel model for R = 2 + 8 with 
8 — 0.1, 0.75, 0.05 and Jo = 1.2. The diamonds indicate the 
results from numerical integration for the 8 = 0.1 case. 



model and the density stratification being determined by: 



-gp'/p : 



Jo 



cosh 2R (y) 



The advantage of this profile is that there are analytic 
solutions for the Kelvin-Helmholtz stability boundaries 
Jo(fc) for the cases that R = 0, 1 and 2. In addition 
there is an analytic solution Jo(k) for the modes k for 
which c(k) = 1 for the R = 2 case. The R = case was 
examined by |4l| where it was shown that the Kelvin 
Helmholtz unstable modes satisfy Jo < fc 2 (l — A: 2 ). The 
R = 1 case (that reduces to the R — 1 case of the Hazel 
model) was investigated by j3^ where it was shown that 
the Kelvin Helmholtz unstable modes satisfy Jo < k(l — 
k). The R = 2 case has not been investigated before 
(to the author's knowledge). One can show follo wing the 
same methods used for the R — 0, 1 cases [3^, El E2 
that the c(k) = modes satisfy: 



Jn = 



k(l - k)(2 + k)(3 + k) 
4(/c + l) 2 



with 4>{y) = [1 - tan%) a ] fc/2 • [tanh(j/)] 1 / 4 -V 1 /4-'o 
and provide the Kelvin-Helmholtz instability boundary. 
c(k) = 1 modes on the other hand that are of interest for 
the Holmboe instability satisfy 



Jn = 



k(3 + 2k) 
(fc + l) 2 



for k < 1. The stream- function 
given by <j> = [1 + tanh(j/)] fc/2 • [1 - 



b for these modes is 
tanh(y)]\ /F ^ 



4+1 -Jo 



The Kelvin-Helmholtz stability boundaries for the three 
cases R = 0, 1, 2 along with the c = 1 solutions for the 
R = 2 case are shown in figure For this example we 
also have that a = 2 and g = fc/2 but J x = 2 2R ~ 4 Jq. 
The Jo(fc) relation for the c = 1 solutions does not sat- 
isfy the criterion[2]that now reads Jo < k/(k + l) 2 , thus 
the R = 2 case is stable (to the examined order) and 
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is the critical value above which the Holmboe instabil- 
ity begins. Because Joo is four times bigger than in the 
Hazel model (for the same Jo) the resulting growth rate 
is smaller by a factor of 4r 2X f s that is close to 10~ 14 for 
the 5 = 0.1 case. Figure El shows the growth rate based 
on the asymptotic expansion 1251 No numerical results 
could be obtained for this case for values of S smaller 
than S < 0.1 that would justify a comparison with the 
asymptotic expansion. This example, when compared 
with one of Hazel, clearly demonstrates the sensitivity of 
the resulting growth rate to the large y asymptotic be- 
havior of J(y) and U(y): a change by a factor of 4 in 
Joo resulted in 14 orders of magnitude difference in the 
growth rate. 



V. CONCLUSIONS 

In this paper we have examined analytically Holmboe's 
instability for smooth density and velocity profiles. We 



8 



have shown for a large family of flows that the modes with 
phase velocity equal to the maximum or minimum of the 
unperturbed velocity profile when they exist and if the 
parameter R is above the critical value (R cr it = 2) they 
constitute a stability boundary. This result confirms the 
results obtained numerically in where the fact that 



the 



max / min 



modes are marginally unstable was 



only conjectured based on physical arguments. It is also 
the first time shown analytically that the value of R = 2 
for the Hazel model is the critical value R cr it above which 
the Holmboe instability begins. 

For the case that the parameter R is only slightly larger 
than it's critical value R cr it = 2, the dispersion relation 
c(k) was obtained based on an asymptotic expansion. 
For this marginally unstable flow the growth rate £ as 
well as the width of the instability stripe Aq have a very 
strong dependence on the deviation of R from it's criti- 
cal value. In particular the growth rate £ and the width 
of the instability Ak = scale as exp[— 2Xj/(R — R cr it)] 
and exp[— 7/(i?— Rcrit)] respectively (for some positive 
constant 7). For this reason the numerical investigations 
performed in the past 0, 0] were not able to capture 
the instability for values of R very close to R cr it- 

We believe also that the present results go beyond the 
clarifying of a mathematical detail in the literature. They 
clearly demonstrate the nature of the Holmboe instabil- 
ity in a quantitative way. The critical layer for Holmboe 
unstable modes appears at large heights where the shear 
rate can overcome stratification and the behavior mode 
will strongly depend on the properties of the shear at this 
height. The critical layer is then coupled to the grav- 
ity wave modes that their properties are determined by 
y = 0(1) quantities and are fast enough to travel with the 
same velocity as the velocity at the critical layer height. 
It is this coupling that gives rise to the instability, and 
it is restricted only to the wave numbers that result in a 
phase speed that is smaller than the maximum velocity 
of the shear but big enough so that the critical layer is at 
a large enough height so that Ri(y) < 1/4. Therefore in 
any experimental setup, precise measurements of the ve- 
locity and density stratification are needed in both large 
and small heights in order to make comparisons of the 
measured growth rate and the theoretical predictions. 

Finally we believe that the results given in this paper 
can provide a basis for further numerical and analytical 
investigations such as an examination of the weakly non- 
linear theory where a small non-linearity is taken into 
account in order to examine the long time evolution of 
an unstable mode beyond the linear stage. 
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APPENDIX A: THE HYPER-GEOMETRIC 
EQUATION: BASIC PROPERTIES 



The Hypergeometric equation is: 

d 2 f df 
z(l - z)-4 + [d - (a + b + l)z]4- - abf = 0. (Al) 
dz A 



dz 



The solution that remains finite as z — > is the Hyperge- 
ometric function: / = F(a, b, c; z). For the normalization 
condition we are using we have the following limits: 

lim F(a, b, c; z) = 1/T(d) 



T(d)T(d-a-b) 
hm F(a,b,c;z) * T{d _ a)r{d _ b) 

T{d)T(a + b- d) 



r(o)r(6) 



(l-z) 



d—a—b 



lim F(a,b,c;z) * W^,- z) -* 



T(d)T(a - b) 
T(a)T(d-b) 



(-*) 



-b 



provided that d 7^ 0, —1, —2, . . . and a — b is not an inte- 
ger. 



APPENDIX B: THE SIGN OF THE INSTABILITY 
TERM 

To determine whether we have successful matching or 
not we need to find the sign of the imaginary part in the 
dispersion relational We examine each term separately. 
Clearly T(a), T(2A) and T(d — b) are all positive factors 
since the argument of the T-function is positive. The fac- 
tor r(— 2A) is changing sign every time 2A is an integer. 
However its product with sin(2A7r) is always remains neg- 
ative. The factors T(b) and T(d~ a) however can change 
sign depending on value of J. Using the expressions for 
a, b, done can show that -1 < b < if J < (2q)/(2q + l) 2 
or if 2q < 1, and positive otherwise. Similarly we have 
that -1 < d-a < if J > (2q)/(2q+l) 2 and 2q < 1 and 
non-negative otherwise. Combining these two inequali- 
ties we can determine the sign of the product 

r(6)r(d-o) < if and only if J < (2q)/(2q + I) 2 . 

The result in 1201 then follows. 

To find the dispersion relation for the small 6 and 
e given by 1221 we need to find an expression for the 
term T(b)T(d — a). Substituting the choice of e given 
by 1221 in the expression for b and d — a and using 
J = Joq{£ c i/v) 5 — ^oo(e) <5 (l + ^l n (ci/c)) w e have that 
to first order in 5 if 2go > 1' 



b ~ -SJoot 5 ln(ci/er) 



1 



1 



VV4- Jooe S yjl + q - Joo£ S 
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and b = 0(1) if 2q < 1. Similarly, 



where 



d—a~ -5Joc€ & ln(ci/er) 



Ci 



2In 



Joot\J {2qoh\ 2X sm(2X)T(a)T(d-b)T(-2X) 



r(to)r(2A) 



if 2g < 1 and d—a = 0(1) if 2q > 1. Using the T- function 
property r(5) = T(l + S)/S we can write the dispersion 
relation for g < q < qo + eah/koh as: 



c = l-2(g-g )go/i//2+<5Ci (qo-q)^ In 



2(g - go)go/i 
Ii<J 



yi/4- Jooe 5 v 71 + % - -W 5 
with tu = 6 if 2g < 1 and to = e? — a if 2g > 1. 
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